Pan-Cancer Analysis of Patient Tumor Single-Cell Transcriptomes Identifies Promising Selective and Safe Chimeric Antigen Receptor Targets in Head and Neck Cancer.

Chimeric antigen receptor (CAR) T cell therapies have yielded transformative clinical successes for patients with blood tumors, but their full potential remains to be unleashed against solid tumors. One challenge is finding selective targets, which we define intuitively to be cell surface proteins that are expressed widely by cancer cells but minimally by healthy cells in the tumor microenvironment and other normal tissues. Analyzing patient tumor single-cell transcriptomics data, we first defined and quantified selectivity and safety scores of existing CAR targets for indications in which they are in clinical trials or approved. We then sought new candidate cell surface CAR targets that have better selectivity and safety scores than those currently being tested. Remarkably, in almost all cancer types, we could not find such better targets, testifying to the near optimality of the current target space. However, in human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSC), for which there is currently a dearth of existing CAR targets, we identified a total of twenty candidate novel CAR targets, five of which have both superior selectivity and safety scores. These newly identified cell surface targets lay a basis for future investigations that may lead to better CAR treatments in HNSC.

Optimizing cancer immunotherapy response prediction by tumor aneuploidy score and fraction of copy number alterations.

Identifying patients that are likely to respond to cancer immunotherapy is an important, yet highly challenging clinical need. Using 3139 patients across 17 different cancer types, we comprehensively studied the ability of two common copy-number alteration (CNA) scores-the tumor aneuploidy score (AS) and the fraction of genome single nucleotide polymorphism encompassed by copy-number alterations (FGA)-to predict survival following immunotherapy in both pan-cancer and individual cancer types. First, we show that choice of cutoff during CNA calling significantly influences the predictive power of AS and FGA for patient survival following immunotherapy. Remarkably, by using proper cutoff during CNA calling, AS and FGA can predict pan-cancer survival following immunotherapy for both high-TMB and low-TMB patients. However, at the individual cancer level, our data suggest that the use of AS and FGA for predicting immunotherapy response is currently limited to only a few cancer types. Therefore, larger sample sizes are needed to evaluate the clinical utility of these measures for patient stratification in other cancer types. Finally, we propose a simple, non-parameterized, elbow-point-based method to help determine the cutoff used for calling CNAs.

Predicting tumor immune microenvironment and checkpoint therapy response of head & neck cancer patients from blood immune single-cell transcriptomics.

The immune state of tumor microenvironment is crucial for determining immunotherapy response but is not readily accessible. Here we investigate if we can infer the tumor immune state from the blood and further predict immunotherapy response. First, we analyze a dataset of head and neck squamous cell carcinoma (HNSCC) patients with matched scRNA-Seq of peripheral blood mononuclear cells (PBMCs) and tumor tissues. We find that the tumor immune cell fractions of different immune cell types and many of the genes they express can be inferred from the matched PBMC scRNA-Seq. Second, analyzing another HNSCC dataset with PBMC scRNA-Seq and immunotherapy response, we find that the inferred ratio between tumor memory B and regulatory T cell fractions is predictive of immunotherapy response and is superior to the well-established cytolytic and exhausted T-cell signatures. Overall, these results showcase the potential of scRNA-Seq liquid biopsies in cancer immunotherapy, calling for their larger-scale testing. Significance: This head and neck cancer study demonstrates the potential of using blood single-cell transcriptomics to (1) infer the tumor immune status and (2) predict immunotherapy response from the tumor immune status inferred from blood. These results showcase the potential of single-cell transcriptomics liquid biopsies for further advancing personalized cancer immunotherapy.

Common evolutionary trajectory of short life-cycle in Brassicaceae ruderal weeds.

Weed species are detrimental to crop yield. An understanding of how weeds originate and adapt to field environments is needed for successful crop management and reduction of herbicide use. Although early flowering is one of the weed trait syndromes that enable ruderal weeds to overcome frequent disturbances, the underlying genetic basis is poorly understood. Here, we establish Cardamine occulta as a model to study weed ruderality. By genome assembly and QTL mapping, we identify impairment of the vernalization response regulator gene FLC and a subsequent dominant mutation in the blue-light receptor gene CRY2 as genetic drivers for the establishment of short life cycle in ruderal weeds. Population genomics study further suggests that the mutations in these two genes enable individuals to overcome human disturbances through early deposition of seeds into the soil seed bank and quickly dominate local populations, thereby facilitating their spread in East China. Notably, functionally equivalent dominant mutations in CRY2 are shared by another weed species, Rorippa palustris, suggesting a common evolutionary trajectory of early flowering in ruderal weeds in Brassicaceae.

Defining the scope for altering rice leaf anatomy to improve photosynthesis: a modelling approach.

Leaf structure plays an important role in photosynthesis. However, the causal relationship and the quantitative importance of any single structural parameter to the overall photosynthetic performance of a leaf remains open to debate. In this paper, we report on a mechanistic model, eLeaf, which successfully captures rice leaf photosynthetic performance under varying environmental conditions of light and CO2 . We developed a 3D reaction-diffusion model for leaf photosynthesis parameterised using a range of imaging data and biochemical measurements from plants grown under ambient and elevated CO2 and then interrogated the model to quantify the importance of these elements. The model successfully captured leaf-level photosynthetic performance in rice. Photosynthetic metabolism underpinned the majority of the increased carbon assimilation rate observed under elevated CO2 levels, with a range of structural elements making positive and negative contributions. Mesophyll porosity could be varied without any major outcome on photosynthetic performance, providing a theoretical underpinning for experimental data. eLeaf allows quantitative analysis of the influence of morphological and biochemical properties on leaf photosynthesis. The analysis highlights a degree of leaf structural plasticity with respect to photosynthesis of significance in the context of attempts to improve crop photosynthesis.

3dCAP-Wheat: An Open-Source Comprehensive Computational Framework Precisely Quantifies Wheat Foliar, Nonfoliar, and Canopy Photosynthesis.

Canopy photosynthesis is the sum of photosynthesis of all above-ground photosynthetic tissues. Quantitative roles of nonfoliar tissues in canopy photosynthesis remain elusive due to methodology limitations. Here, we develop the first complete canopy photosynthesis model incorporating all above-ground photosynthetic tissues and validate this model on wheat with state-of-the-art gas exchange measurement facilities. The new model precisely predicts wheat canopy gas exchange rates at different growth stages, weather conditions, and canopy architectural perturbations. Using the model, we systematically study (1) the contribution of both foliar and nonfoliar tissues to wheat canopy photosynthesis and (2) the responses of wheat canopy photosynthesis to plant physiological and architectural changes. We found that (1) at tillering, heading, and milking stages, nonfoliar tissues can contribute ~4, ~32, and ~50% of daily gross canopy photosynthesis (A cgross; ~2, ~15, and ~-13% of daily net canopy photosynthesis, A cnet) and absorb ~6, ~42, and ~60% of total light, respectively; (2) under favorable condition, increasing spike photosynthetic activity, rather than enlarging spike size or awn size, can enhance canopy photosynthesis; (3) covariation in tissue respiratory rate and photosynthetic rate may be a major factor responsible for less than expected increase in daily A cnet; and (4) in general, erect leaves, lower spike position, shorter plant height, and proper plant densities can benefit daily A cnet. Overall, the model, together with the facilities for quantifying plant architecture and tissue gas exchange, provides an integrated platform to study canopy photosynthesis and support rational design of photosynthetically efficient wheat crops.

Diurnal and Seasonal Variations of Photosynthetic Energy Conversion Efficiency of Field Grown Wheat.

Improving canopy photosynthetic light use efficiency and energy conversion efficiency (ε c ) is a major option to increase crop yield potential. However, so far, the diurnal and seasonal variations of canopy light use efficiency (LUE) and ε c are largely unknown due to the lack of an efficient method to estimate ε c in a high temporal resolution. Here we quantified the dynamic changes of crop canopy LUE and ε c during a day and a growing season with the canopy gas exchange method. A response curve of whole-plant carbon dioxide (CO2) flux to incident photosynthetically active radiation (PAR) was further used to calculate ε c and LUE at a high temporal resolution. Results show that the LUE of two wheat cultivars with different canopy architectures at five stages varies between 0.01 to about 0.05 mol CO2 mol-1 photon, with the LUE being higher under medium PAR. Throughout the growing season, the ε c varies from 0.5 to 3.7% (11-80% of the maximal ε c for C3 plants) with incident PAR identified as a major factor controlling variation of ε c . The estimated average ε c from tillering to grain filling stages was about 2.17%, i.e., 47.2% of the theoretical maximal. The estimated season-averaged radiation use efficiency (RUE) was 1.5-1.7 g MJ-1, which was similar to the estimated RUE based on biomass harvesting. The large variations of LUE and ε c imply a great opportunity to improve canopy photosynthesis for greater wheat biomass and yield potential.

Potential metabolic mechanisms for inhibited chloroplast nitrogen assimilation under high CO2.

Improving photosynthesis is considered a major and feasible option to dramatically increase crop yield potential. Increased atmospheric CO2 concentration often stimulates both photosynthesis and crop yield, but decreases protein content in the main C3 cereal crops. This decreased protein content in crops constrains the benefits of elevated CO2 on crop yield and affects their nutritional value for humans. To support studies of photosynthetic nitrogen assimilation and its complex interaction with photosynthetic carbon metabolism for crop improvement, we developed a dynamic systems model of plant primary metabolism, which includes the Calvin-Benson cycle, the photorespiration pathway, starch synthesis, glycolysis-gluconeogenesis, the tricarboxylic acid cycle, and chloroplastic nitrogen assimilation. This model successfully captures responses of net photosynthetic CO2 uptake rate (A), respiration rate, and nitrogen assimilation rate to different irradiance and CO2 levels. We then used this model to predict inhibition of nitrogen assimilation under elevated CO2. The potential mechanisms underlying inhibited nitrogen assimilation under elevated CO2 were further explored with this model. Simulations suggest that enhancing the supply of α-ketoglutarate is a potential strategy to maintain high rates of nitrogen assimilation under elevated CO2. This model can be used as a heuristic tool to support research on interactions between photosynthesis, respiration, and nitrogen assimilation. It also provides a basic framework to support the design and engineering of C3 plant primary metabolism for enhanced photosynthetic efficiency and nitrogen assimilation in the coming high-CO2 world.

Scheduling of nitrogen fertilizer topdressing during panicle differentiation to improve grain yield of rice with a long growth duration.

Topdressing at panicle differentiation (PF) according to soil fertility and regularity of rice nutrient absorption is an important agronomic practice used in cultivation of rice cultivars with a long growth duration. We studied the impacts of timing of nitrogen fertilizer application during PF on photosynthesis and yield-related agronomic traits in 'Y-Liang-You 900' and 'Y-Liang-You 6', which are representative rice cultivars with a long growth duration. Data for two years showed that timing of topdressing application during PF affected panicles per unit area, percentage grain set, spikelets per panicle, and leaf photosynthetic traits during the grain-filling period. Topdressing at the initial stage of flag-leaf extension resulted in higher grain yield (typically by 10.55-19.95%) than in plants without topdressing. Grain yield was significantly correlated with flag leaf photosynthetic rate and leaf SPAD value (r = 0.5640 and r = 0.5589, respectively; p < 0.01) at an advanced grain-filling stage (30 days after heading). Surprisingly, grain yield was not correlated with carbohydrate remobilization from the stem and sheath. For rice cultivars with a long growth duration, nitrogen-fertilizer topdressing must be applied at the initial stage of flag-leaf extension to delay leaf senescence during the grain-filling stage and realize the enhanced yield potential.

Morphological and physiological factors contributing to early vigor in the elite rice cultivar 9,311.

Huanghuazhan (HHZ) and 9,311 are two elite rice cultivars in China. They have achieved high yield through quite different mechanisms. One of the major features that gives high yield capacity to 9,311 is its strong early vigor, i.e., faster establishment of its seedling as well as its better growth in its early stages. To understand the mechanistic basis of early vigor in 9,311, as compared to HHZ the cultivar, we have examined, under controlled environmental conditions, different morphological and physiological traits that may contribute to its early vigor. Our results show that the fresh weight of the seeds, at germination, not only determined the seedling biomass at 10 days after germination (DAG), but was also responsible for ~ 80% of variations in plant biomass between the two cultivars even up to 30 DAG. Furthermore, the 9,311 cultivar had a larger root system, which led to its higher nitrogen uptake capacity. Other noteworthy observations about 9,311 being a better cultivar than HHZ are: (i) Ten out of 15 genes involved in nitrogen metabolism were much more highly expressed in its roots; (ii) it had a higher water uptake rate, promoting better root-to-shoot nitrogen transfer; and (iii) consistent with the above, it had higher leaf photosynthetic rate and stomatal conductance. All of the above identified features explain, to a large extent, why the 9,311, as compared to HHZ, exhibits much more vigorous early growth.

Architectural and Physiological Features to Gain High Yield in an Elite Rice Line YLY1.

Identification of traits strongly associated with high yield can help future gene engineering towards improvements of productivity. Here we systematically determine the major architectural and physiological features associated with high yield in two elite historical hybrid rice cultivars, i.e., YLY1 and LYP9. Data from a six-year experiment show that high yield of YLY1 are related to a number of architectural and physiological parameters. Compared to LYP9, YLY1 had 5.5% and 47.3% higher canopy photosynthesis under high and low photosynthetic photon flux densities, respectively, during the grain filling stage, an average 1.5% higher proportion of biomass allocation to above-ground tissues, a 4.5%-10.5% higher photosynthate reserve in leaf sheath before grain filling, and a more efficient photosynthate translocation during grain filling and finally an average 25.2% higher number of productive tillers. These features differ dramatically from features associated with high yield in YLY900 and Yongyou12#, two other high-yielding rice cultivars in China. These identified features and their combinations can support designing new strategies in the future high-yield rice breeding.

An in situ approach to characterizing photosynthetic gas exchange of rice panicle.

BACKGROUND: Photosynthesis of reproductive organs in C3 cereals is generally regarded as important to crop yield. Whereas, photosynthetic characteristics of reproductive organs are much less understood as compared to leaf photosynthesis, mainly due to methodological limitations. To date, many indirect methods have been developed to study photosynthesis of reproductive organs and its contribution to grain yield, such as organ shading, application of herbicides and photosynthetic measurement of excised organs or tissues, which might be intrusive and cause biases. Thus, a robust and in situ approach needs to be developed. RESULTS: Here we report the development of a custom-built panicle photosynthesis chamber (P-chamber), which can be connected to standard infrared gas analyzers to study photosynthetic/respiratory rate of a rice panicle. With the P-chamber, we measured panicle photosynthetic characteristics of seven high-yielding elite japonica, japonica-indica hybrid and indica rice cultivars. Results show that, (1) rice panicle is photosynthetically active during grain filling, and there are substantial inter-cultivar variations in panicle photosynthetic and respiratory rates, no matter on a whole panicle basis, on an area basis or on a single spikelet basis; (2) among the seven testing cultivars, whole-panicle gross photosynthetic rates are 17-54 nmol s-1 5 days after heading under photon flux density (PFD) of 2000 µmol (photons) m-2 s-1, which represent some 20-38% of that of the corresponding flag leaves; (3) rice panicle photosynthesis has higher apparent CO2 compensation point, light compensation point and apparent CO2 saturation point, as compared to that of a typical leaf; (4) there is a strong and significant positive correlation between gross photosynthetic rate 5 days after heading on a single spikelet basis and grain setting rate at harvest (Pearson correlation coefficient r = 0.93, p value < 0.0001). CONCLUSIONS: Rice panicle gross photosynthesis is significant, has great natural variation, and plays an underappreciated role in grain yield formation. The P-Chamber can be used as a tool to study in situ photosynthetic characteristics of irregular non-foliar plant organs, such as ears, culms, leaf sheaths, fruits and branches, which is a relatively less explored area in current cereal breeding community.

Genome-Wide Association Study Unravels LRK1 as a Dark Respiration Regulator in Rice (Oryza sativa L.).

Respiration is a major plant physiological process that generates adenosine triphosphate (ATP) to support the various pathways involved in the plant growth and development. After decades of focused research on basic mechanisms of respiration, the processes and major proteins involved in respiration are well elucidated. However, much less is known about the natural variation of respiration. Here we conducted a survey on the natural variation of leaf dark respiration (Rd) in a global rice minicore diversity panel and applied a genome-wide association study (GWAS) in rice (Oryza sativa L.) to determine candidate loci associated with Rd. This rice minicore diversity panel consists of 206 accessions, which were grown under both growth room (GR) and field conditions. We found that Rd shows high single-nucleotide polymorphism (SNP) heritability under GR and it is significantly affected by genotype-environment interactions. Rd also exhibits strong positive correlation to the leaf thickness and chlorophyll content. GWAS results of Rd collected under GR and field show an overlapped genomic region in the chromosome 3 (Chr.3), which contains a lead SNP (3m29440628). There are 12 candidate genes within this region; among them, three genes show significantly higher expression levels in accessions with high Rd. Particularly, we observed that the LRK1 gene, annotated as leucine rich repeat receptor kinase, was up-regulated four times. We further found that a single significantly associated SNPs at the promoter region of LRK1, was strongly correlated with the mean annual temperature of the regions from where minicore accessions were collected. A rice lrk1 mutant shows only ~37% Rd of that of WT and retarded growth following exposure to 35 °C for 30 days, but only 24% reduction in growth was recorded under normal temperature (25 °C). This study demonstrates a substantial natural variation of Rd in rice and that the LRK1 gene can regulate leaf dark respiratory fluxes, especially under high temperature.

Overexpression of maize transcription factor mEmBP-1 increases photosynthesis, biomass, and yield in rice.

Identifying new options to improve photosynthetic capacity is a major approach to improve crop yield potential. Here we report that overexpression of the gene encoding the transcription factor mEmBP-1 led to simultaneously increased expression of many genes in photosynthesis, including genes encoding Chl a,b-binding proteins (Lhca and Lhcb), PSII (PsbR3 and PsbW) and PSI reaction center subunits (PsaK and PsaN), chloroplast ATP synthase subunit, electron transport reaction components (Fd1 and PC), and also major genes in the Calvin-Benson-Bassham cycle, including those encoding Rubisco, glyceraldehyde phosphate dehydrogenase, fructose bisphosphate aldolase, transketolase, and phosphoribulokinase. These increased expression of photosynthesis genes resulted in increased leaf chlorophyll pigment, photosynthetic rate, biomass growth, and grain yield both in the greenhouse and in the field. Using EMSA experiments, we showed that mEmBP-1a protein can directly bind to the promoter region of photosynthesis genes, suggesting that the direct binding of mEmBP-1a to the G-box domain of photosynthetic genes up-regulates expression of these genes. Altogether, our results show that mEmBP-1a is a major regulator of photosynthesis, which can be used to increase rice photosynthesis and yield in the field.

A three-dimensional canopy photosynthesis model in rice with a complete description of the canopy architecture, leaf physiology, and mechanical properties.

In current rice breeding programs, morphological parameters such as plant height, leaf length and width, leaf angle, panicle architecture, and tiller number during the grain filling stage are used as major selection targets. However, so far, there is no robust approach to quantitatively define the optimal combinations of parameters that can lead to increased canopy radiation use efficiency (RUE). Here we report the development of a three-dimensional canopy photosynthesis model (3dCAP), which effectively combines three-dimensional canopy architecture, canopy vertical nitrogen distribution, a ray-tracing algorithm, and a leaf photosynthesis model. Concurrently, we developed an efficient workflow for the parameterization of 3dCAP. 3dCAP predicted daily canopy RUE for different nitrogen treatments of a given rice cultivar under different weather conditions. Using 3dCAP, we explored the influence of three canopy architectural parameters-tiller number, tiller angle and leaf angle-on canopy RUE. Under different weather conditions and different nitrogen treatments, canopy architecture optimized by manipulating these parameters can increase daily net canopy photosynthetic CO2 uptake by 10-52%. Generally, a smaller tiller angle was predicted for most elite rice canopy architectures, especially under scattered light conditions. Results further show that similar canopy RUE can be obtained by multiple different parameter combinations; these combinations share two common features of high light absorption by leaves in the canopy and a high level of coordination between the nitrogen concentration and the light absorbed by each leaf within the canopy. Overall, this new model has potential to be used in rice ideotype design for improved canopy RUE.

Systems models, phenomics and genomics: three pillars for developing high-yielding photosynthetically efficient crops.

Recent years witnessed a stagnation in yield enhancement in major staple crops, which leads plant biologists and breeders to focus on an urgent challenge to dramatically increase crop yield to meet the growing food demand. Systems models have started to show their capacity in guiding crops improvement for greater biomass and grain yield production. Here we argue that systems models, phenomics and genomics combined are three pillars for the future breeding for high-yielding photosynthetically efficient crops (HYPEC). Briefly, systems models can be used to guide identification of breeding targets for a particular cultivar and define optimal physiological and architectural parameters for a particular crop to achieve high yield under defined environments. Phenomics can support collection of architectural, physiological, biochemical and molecular parameters in a high-throughput manner, which can be used to support both model validation and model parameterization. Genomic techniques can be used to accelerate crop breeding by enabling more efficient mapping between genotypic and phenotypic variation, and guide genome engineering or editing for model-designed traits. In this paper, we elaborate on these roles and how they can work synergistically to support future HYPEC breeding.

Systems model-guided rice yield improvements based on genes controlling source, sink, and flow.

A large number of genes related to source, sink, and flow have been identified after decades of research in plant genetics. Unfortunately, these genes have not been effectively utilized in modern crop breeding. This perspective paper aims to examine the reasons behind such a phenomenon and propose a strategy to resolve this situation. Specifically, we first systematically survey the currently cloned genes related to source, sink, and flow; then we discuss three factors hindering effective application of these identified genes, which include the lack of effective methods to identify limiting or critical steps in a signaling network, the misplacement of emphasis on properties, at the leaf, instead of the whole canopy level, and the non-linear complex interaction between source, sink, and flow. Finally, we propose the development of systems models of source, sink and flow, together with a detailed simulation of interactions between them and their surrounding environments, to guide effective use of the identified elements in modern rice breeding. These systems models will contribute directly to the definition of crop ideotype and also identification of critical features and parameters that limit the yield potential in current cultivars.

Source-sink interaction: a century old concept under the light of modern molecular systems biology.

Many approaches to engineer source strength have been proposed to enhance crop yield potential. However, a well-co-ordinated source-sink relationship is required finally to realize the promised increase in crop yield potential in the farmer's field. Source-sink interaction has been intensively studied for decades, and a vast amount of knowledge about the interaction in different crops and under different environments has been accumulated. In this review, we first introduce the basic concepts of source, sink and their interactions, then summarize current understanding of how source and sink can be manipulated through both environmental control and genetic manipulations. We show that the source-sink interaction underlies the diverse responses of crops to the same perturbations and argue that development of a molecular systems model of source-sink interaction is required towards a rational manipulation of the source-sink relationship for increased yield. We finally discuss both bottom-up and top-down routes to develop such a model and emphasize that a community effort is needed for development of this model.